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ABSTRACT 

We have carried out an investigation of the chemical evolution of gas in different carbon- 
rich circumstellar environments. Previous studies have tended to invoke terrestrial flame 
chemistries, based on acetylene (C2H2) combustion to model the formation of carbon dust, via 
Polycyclic Aromatic Hydrocarbons (PAHs). In this work we pay careful attention to the accu- 
rate calculation of the molecular photoreaction rate coefficients to ascertain whether there is 
a universal formation mechanism for carbon dust in strongly irradiated astrophysical environ- 
ments. A large number of possible chemical channels may exist for the formation of PAHs, so 
we have concentrated on the viability of the formation of the smallest building block species, 
C2H2, in a variety of carbon-rich stellar outflows. C2H2 is very sensitive to dissociation by 
UV radiation. This sensitivity is tested, using models of the time-dependent chemistry. We find 
that C2H2 formation is sensitive to some of the physical parameters and that in some known 
sources of dust-formation it can never attain appreciable abundances. Therefore multiple (and 
currently ill-defined) dust-formation channels must exist. 
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1 INTRODUCTION 

The obsei-vational properties of dust and the physical consequences 
of its presence in the inters tellar medium and extremely well- 
known and well-documented ( lDrainel2003l) . However its composi- 
tion, structure and size-distribution are still subjects of much dis- 
cussion. Even though there are no direct samples of interstellar 
dust, carbonaceous dust deposits found in meteorites give us an in- 
dication of what structure and composition we might expect. These 
grains consist of a core of microcrystalline aromatic moi eties sur- 
rounded by a crystalline graphite mantle (lTielenslll997l) . Signifi- 
cantly, the precise chemical kinetic pathways leading to the (rel- 
atively simple) case of carbon-dust formation are almost entirely 
speculative. 

The recent observations of h igher than expected dust masses 
around high redsh ift quasars ( IDwek. Galliano & JonesI l2007t 



ISertoldi et al.ll2003h and local supemovae JMatsuura et al.ll201 If) 
has caused a resurgence of interest in the issues of dust formation 
mechansims and efficiencies. In the local universe, dust is observed 
mainly in two environments: dense interstellar clouds and circum- 
stellar media. The main sources of interstell ar dust are believed to 
be ev olve d red giants, n ovae and supemovae ( lDwe]a2004i ISalpeterl 
1 19771 : see lTielensluOOSl for a review). There are thus two questions 
that we address in this paper: (i) what is the chemical pathway for 
the formation of macromolecular species that leads to the formation 
of dust nucleation sites?, and (ii) what range of physical parameters 
determine the viability of these pathways? 
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molecular processes - stars; AGB and post-AGB - dust, extinc- 



For the sake of (relative) simplicity we limit our discussion 
to the formation of carbon dust in carbon-rich circumstellar envi- 
ronments, making the simplifying assumption that the oxygen is 
completely locked up in the relatively inert form of CO. In such an 
environment the C:0 ratio - 1 is thus a crucial parameter. However, 
it should be noted (see below) that this usually-assumed approxi- 
mation has been challenged by relatively recent observational and 
modelling efforts. 

Some of the earliest attempts at modelling carbon dust for- 
mation concentrated on novae, some of which are seen to pro- 
duce very optically t hick (carbon) dust shells on timescales of 
davs. lRawlingsl jl988r) developed a proto-PDR model of the chem- 
istry in novae ejectae, and found that molecule formation requires 
the presence of H2 and carbon neutrality (CI). The presence of 
C+ (and the carbon ionization continuum) inhibits molecule for- 
mati on. Chemical equilibrium is established very fast (within ~ 
10s). lRawlings & William"sl l ll989l) , hereafter RW89, extended this 
and proposed a non-LTE chemical kinetic model of (carbon) dust 
formation in novae ejectae which proceeds via a limited hydro- 
carbon chemistry - totally dominated by the photophysics, and 
requiring all oxygen to be locked up in CO. This latter find- 
ing has been challenged - both by the observations that suggest 
that carbon dust is formed in oxygen ric h (0>C) environments 
(e.g. lAndrea. Drechsel & Starrfie ld 1994) - and also model s of 
the chemistry in nova outflows (Pontefrac t & Rawlingsll2004) . In 
these models C2 (for which the more significant of the two pre- 
dissociation bands is shielded against photodissociation in a CI re- 
gion) is the key building block species, rather than C2H2 which is 
susceptible to photodissociation. As such, nucleation of dust was 



2 H. Dhanoa & J.M.C. Rawlings 



found to be possible only in a CI region, but requires some source 
of io nization to drive the chemistry . 

IPontefract & Rawlingsl ( l2004h extended the worlc by RW89 
by including oxygen and nitrogen chemistry and updated the pho- 
toreaction rates. Surprisingly, they found that the overall chemistry 
is not as dependent on the photochemistry as earlier studies had 
suggested and that carbon dust could be formed even when all of 
the oxygen is not 'locked' in CO. The work highlighted that CN, 
CO and Ct are key molecular species in the formation pathways of 
larger species. Although C2H2 was included its chemistiy was not 
closely scrutinised. 

More generally, it has been proposed that Polycyclic Aro- 
matic Hydrocarbon molecules (PAHs) are key intermediates in 
carbon dust formation and act as nucleation sites for d ust grains 
faelhng et aLlll996l : ICherchneff. Barker & Tielenslll992h . This in- 
ferred link is due to the the close relationship between the molec- 
ular structure of PAHs and the carbon structure of graphite. Also, 
observations of the far-IR continu um due to car bon dust are ac- 
companied by strong PAH features ( lTielensl2008h . However, some 
caution should be exercised with the interpretation of these fea- 
tures. For example, PAH (and/or hydrogenated amorphous carbon) 
emission is also detected in some carbon-dust producing novae but 
in all cases it is only seen well after dust formation and in an en- 
vironment when it would be susceptible to photodissociation. This 
indicates that - for novae at least - the PAHs are transitory prod- 
ucts of dust des truction, rather than intermediates to dust formation 
JEvans & Rawlings 2008) . 



Gail & Sedlmavn ( ll987D proposed chemical mechanisms re- 



quired to form PAH molecules in stellar winds. The acetylene 
molecule (C2H2) was found to be crucial in the formation path- 
way. A n important contribution came from lFrenklach & FeigelsonI 
( 119891) . hereafter (FF89). They investigated a mechanism of PAH 
formation, using a chemical kinetic approach. The chemistry net- 
work was based on soot production in terrestrial hydrocarbon 
flames, which they applied to astrophysical conditions, i.e. much 
lower pressure and densities. They confirmed that the complex net- 
work of chemical reactions can essentially be described by a se- 
quence of hydrogen abstraction followed by acetylene addition. 

Acetylene does not possess a permanent dipole moment 
and so does not have a rotational spectrum observable at radio 
wavelemgths. It is therefore somewhat elusive; - detections rely 
on rovibrational transitions, observable in the mid-infrared. For 
this reason, detections tend to be limited to the warm gas (^ 
100 — lO OOK) along the lines of sight towards young stellar ob- 
jects (e.g. lBoonman et al.ll2003r) . There is also strong evidence for 
the presence of C2H2 in evolved carbon stars; the most well known 
of which - IR C-I-102I6 - exhibits a forest of C2H2 lines in the 1 1- 
14^m window JFonfn'a et al.l2008h . Indeed, after H2 and CO, C2H2 
is determined to be the most abundant species in the gas - with 
a fractional abundance of ~ 8 x 10^^ - although this is primarily 
detected outside the dust-formation zone and so is not necessar- 
ily associated with the dust-formation process. There may even 
be some evidence of C2H2 freeze-out into icy mantles at larger 
radii in IRC-l-10216, and an icy mantle-based origin for interstellar 
C2H2 is strongly suggested by Spitzer observations of the corre- 
lation between C2H2 and gas-phase CO2 towards C epheus A East 
JSon nentrucker, Gonzalez- Alfonso & Ne ufeldluOOVh . 

Cherchneff, Barker & Tielens ( 1992) addressed the issue of 
PAH chemistry and dust formation in the carbon-rich envelopes 
of late AGB stars; driven by stellar pulstations and (strong, den- 
sity enhancing) shocks. Although extremely dense, the gas is much 
cooler and not subject to the intense FUV radiation field that is 



present in a nova, so the chemical model was based on a terres- 
trial neutral gas acetylenic-buming soot chemistry. In such envi- 
ronments the carbon budget is largely split between CO and C2H2. 
The chemistry was largely based on well-studied flame chemistries 
(eg. FF89), expanded and augmented to include free radicals. It 
was found that C2H2 is the key species in macro-molecular growth 
and (especially) ring closure/cyclization to form aromatic hydro- 
carbon rings, as well as the polymerization to multiple ring PAHs. 
The formation of the first aromatic ring is often recognised as the 
main 'bottleneck' in macro -molecular growth/dust n ucleation. In 
this, as in later studies (e.g. ICherchneff & Caulll999l) . ring closure 
was found to proceed following reactions of two propargyl radi- 
cals (C3H3) or via a lesser channel involving reaction of acetylene 
with 1,3 butadiyne (C4H2) (e.g. see Figs. 3 & 4 in that paper). In 
either case C2H2 initiates the chemistry. Reaction networks need 
to include both two-body and three-body reactions, but did not in- 
clude photoreactions, or an ion-molecule chemistry. Although they 
were able to identify the presence of a temperature window (900- 
llOOK) in which the formation of PAHs is possible, the formation 
efficiency is too low - possibly due to the lack of inclusion of the 
effects of a local UV radiation field and/or an ion-molecule chem- 
istry. 

Considering a different class of object IWoods et al.l ( l2002h 
modelled the formation of benzene (CgHg) in a protoplanetary neb- 
ula (CRL618) where the chemistry is characterised by high den- 
sities, temperatures and ionization rates - although photoreactions 
were significantly inhibited throughout the molecule-formation 
zone due to the assumed presence of dust. This chemistry is some- 
what different to the Cherchneff et al. schemes due to the pres- 
ence of a degree of ionization and an efficient ion-molecule chem- 
istry. However, even in these conditions, C2H2 plays a crucial role 
- mainly through the initiating reaction of 

HCO+ +C2H2 -> products... 

All of the above models and reaction schemes have made extremely 
simplistic assumptions about the radiation field and the photoreac- 
tion rates; either by only considering the chemistry in low luminos- 
ity environments (or external to an optically thick dust-formation 
zone) - in which case they are ignored, or by simply 'scaling up' 
interstellar photoreaction rates in a way that does not properly take 
into account the specific nature of the radiation field and/or the sig- 
nificance of ionization continua. The overall aim of this work is 
therefore to include an ion-molecule chemistry and a careful re- 
calculation of the photoreaction rates using local radiation fields 
and accurate cross-sectional data to re-assess the viability of C2H2 
formation in a variety of carbon-rich circumstellar environments. 
As C2H2 is fundamental to the formation of PAH molecules, we 
try to determine if the amount of C2H2 produced in carbon rich 
AGB and nova environments is enough to seed dust formation. We 
use a specially adapted chemical network for high temperatures and 
densities, paying particular attention to a more accurate account of 
the photochemistry involved in non-interstellar conditions, which 
(as explained above) previously had been highly over-simplified. 
We investigate a broad range of dust-forming objects which in- 
cludes carbon stars with effective temperatures higher than AGB 
stars (e.g. CH subgiants), novae and R Coronae Borealis stars. By 
doing this we aim to establish the range of parameter space within 
which effective C2H2 formation (and by inference PAH foiTnation) 
is possible. 

In section ^ we describe the model that was implemented. 
Section ^specifies the physical and chemical parameters used for 
the three different carbon rich circumstellar environments that we 
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have modelled. We describe the results and analysis in section ^ 
Finally, in section ^we summarise our results and their implica- 
tions. 



2 THE MODEL 

The nucleation sites for dust grain formation are formed in the gas 
phase in a coolin g, expanding atmosp here, such as in the outflow 
from a red giant JCadwell et al.lll994f) . Therefore if PAHs are the 
key intermediate molecules to carbonaceous dust they should also 
be formed in the stellar winds that feed the circumstellar environ- 
ment. We therefore model the time-dependent chemistry that oc- 
curs in an appropriately located gas parcel. This will be the part 
of the envelope whose temperature lies within the range that is ap- 
propriate for dust formation - corresponding to the 'condensatio n 
radius' in previous studies (eg. lClavton & Wickramasinghelll976r) . 

As explained above, the main chemical pathway to form 
a PAH molecule sta rts with the polymerization of C2H2 
JGail & Sedlmavril 19871) and we therefore identify acetylene as be- 
ing the key molecule in the formation of larger species, such as 
PAHs. We investigate the viability of acetylene formation as being 
the controlling factor in the formation of dust nucleation sites via 
PAHs as intermediate macro-molecules. 

The gas in stellar outflows is dense (n> 10 cm^) and 
warm (T> lOOOK), so that (providing the gas is not too hot, 
such that bond-breaking, collisional dissociation reactions can oc- 
cur), molecules can be formed efficiently and quickly. As the 
gas expands and its density falls, the chemical reaction rates are 
quenched. However, the geometrical dilution timescale is typi- 
cally many orders of magnitude larger than the chemical kinetic 
timescale (Rawlings 1988) and this is not a critical factor in deter- 
mining the dust-formation efficiency. Despite the possible presence 
of shocks etc. - and although the gas is very far from thermal equi- 
librium, the characteristic timescale of the outflow is much larger 
than the chemical kinetic timescale. We therefore assume that both 
the circumstellar environment and the stellar wind are spherically 
symmetric around the star and that the physical conditions (density, 
temperature, radiation dilution factor etc.) are not time dependent 
and remain constant. 

The chemistry is evolved to chemical equilibrium which, as 
argued above, occurs on a timescale that is very much shorter than 
that for which changes in the physical conditions occur. In practice, 
equilibrium was determined from the requirement that the abun- 
dances of all species changed by < 1% on any time step and typi- 
cally was achieved on timescales of seconds. We have investigated 
the sensitivity of the equilibrium abundances to the physical con- 
ditions. The key free parameters that we have varied are: the total 
abundances of the elements, the photospheric temperature of the 
star, the density and the temperature of the gas. Other parameters, 
such as the dilution factor for the radiation field, are taken to be 
defined by the observed location of dust formation. The elemental 
abundances are given in Table 4. The cosmic ray ionization rate 
is set to 10^ s^\ which is representative of the assumed galactic 
background rate (iDalgamoii2006i) . 



2.1 The Chemistry 

The chemical network includes 1537 reactions between 102 gas- 
phase species involving 8 elements (H, He, C, N, O, S, Na and Si). 
These are listed in Table 1 . 

The chemical complexity of the reactants was limited as we 



Table 1. Chemical Species 



c, c+,c- 



H, H2 
, CO, CO+, 



CH, CH+, CH2, CHj, CH3, CH+ 



CH4, CH|, C2, Cj, C2H, C2H+, C2H2, C2H^ 



, NH2,NHJ,NH3,NH+, 



NH+ N2, N+ N2H 



2 ' 
H3O+ 



N, N+,NH, NH+ 

O, 0+, O, O2, Oj, OH, OH, 0H+, H2O, H2O 

HCO, HCO+, H2CO, H2CO+, CO2, CO+ 

CN, CN+, CN", HCN, HCN+, HNC, HCNH+ 

NO, N0+, HNO, HNO+, NO2, NO|, OCN 

He, He+, HeH+, Na, Na+, e" 

S, S+, S", HS, HS+, H2S, H2S+, H3S+ 

CS, CS+, HCS, HCS+, H2CS+, SO, SO+, SO2, SOj, OCS, OCS^ 

NS, Si, Si+, SiH, SiH+, SiH+ SiO, SiO+, SiOH+, Si02 



Table 2. Main formation reactions of C2H2, rates from IWoodall et al.l 



i2007l) 



Reaction 



Reaction rate 



1 


CH2 + CH2 -^ C2H2 +2H 


ki=1.80x \Q-^o^-mKlT 


2 


CH2 + CH2 ^ C2H2 +H2 


k2 = 2.63xl0-''e-«i3'^/^ 


3 


H- + C2H ^ C2H2 + e- 


k3 = 1.00xlO-' 


4 


C- + CH2 -^ C2H2 + e- 


k(=5.00x 10-'" 


5 


C + CH3 ^ C2H2 + e- 


kg = 1.00x10-1" 



only concentrate on the initiating steps in the formation of larger 
molecular species. 

The ratefile is an adapted set of reactions, applicable to 
high densities, high temperatures and intense radiation fields. 
Many of the reactions and rate coefficients (k, ) were taken 
from the UMIST database for Astrochemistry (IWoodall et al.l 
120071) where the formulation for the rate coefficients is valid 
in the operative temperature range. Additional data was taken 
from databases that are re levant to h igh tem peratu r e/high 
density environments (e.g. iRawlings. Dr ew & B arlowl Il993t 
Lim. Raga. Rawlings & Williams! 120021 : IPontefract & Rawlingj 
20041 etc.). Although the densities are very much larger than those 
applicable to interstellar studies, they are not so high (Js lO'^ 
cm- ) that three-body reactions need to be taken into account. 

Reaction types that are included in the chemistry are: charge 
transfer, ion-molecule, radiative association, neutral exchanges, ra- 
diative recombination, dissociative recombination and negative ion 
reactions. Reactions which are particularly significant in addition to 
those applicable to interstellar clouds are neutral-neutral reactions, 
reactions with significant activation energy barriers, and reactions 
between ro-vibrationally excited species, due to the high tempera- 
tures that exist in circumstellar outflows. 

The main formation and destruction reactions for C2H2 are 
given in Tables 2 and 3. 



2.1.1 Photoreactions 

In circumstellar environments, irradiated by the central star, the UV 
radiation field is typically very much stronger and has quite a differ- 
ent spectral shape to that whch pertains in the interstellar medium. 
It is therefore important that we pay particular attention to the cor- 
rect evaluation of the photo-reaction rate coefficients: The rates for 
all photodissociation, photoionisation and photodetachment reac- 
tions are calculated ab initio. 

In these conditions, the key factor in determining the viabil- 
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Table 3. Main destruction reactions of C2H2, rates from lWoodall et al.1 1 120071) 



Reaction 



Reaction rate 



leg = 1.45x10- '0(4) ' 
k7 = 1.45xl0-'0(4)-«-'^ 
kg = 5.40 X 10-'" 



9 
10 

11 
12 
13 
14 
15 
16 
17 
18 
19 
20 



C + C2H2 
C + C2H2 
H+ C2H2 - 
H + C2H2 
O + C2H2 
O + C2H2 
He+ +C2H2 -^ 
He+ +C2H2 - 
He+ +C2H2 ^ 
He+ +C2H2 
C+ + C2H2 
S+ + C2H2 
Si+ + C2H2 
Si+ + C2H2 - 
Si + C2H 



► C3H + H 

> C3 + H2 

^ C2H+ + H 

-^ C2H +H2 
^ CO + CH2 
-^ C2H +0H 
. CH+ + CH +He 

-^Ct + He + H2 
• C2H+ + He + H 
-> C2HJ + He 
^C3H+ + H 
-^ HC2S+ H 
-^ SiC2H+ + H 
SiC2Hj + photon 
-^ SICt +H7 



-0.12 



k9 = 3.80xl0-"'e-^ 



kio = 8.39xlO-'2(3lf5) 
kii =5.30x lO-'e- 
ki2=7.70x 10-'" 
ki3 = 1.61 X K)-' 
ki4 = 8.75x 10-'° 
ki5=2.54x 10-'° 
ki6 = 2.20x 10-5 
ki7=9.50x 10-'° 

ki8 = 1.30xlO-'°(5l5)-''" 
ki9=2.50x 10-'° 

k20=1.00xl0-'3(3^; 



1.03 1197 

e-^7- 

8520 



29 



1.00 



ity of any chemical kinetic pathway that leads to dust-formation is 
therefore the stability of the intermediates against photodissocia- 
tion/photoionisation by the stellar radiation field. 

The calculation of the photorates by a time-dependent non- 
interstellar radiation field in an expanding atmosphere is a non- 
trivial task. There are two components to the calculation of the 
photoreaction rate-coefficients in the local radiation field: 

(i) a basic 'photospheric' rate, calculated for the specific pho- 
toreaction and radiation field 

(ii) a time-dependent geometrical dilution, due to the expansion 
of the wind 

In earlier models, and for the sake of simplicity, photoreac- 
tion rates were calculated by scaling the unshielded interstellar 
photorates by a ratio of the stellar flux to that of the interstellar 
radiation field, either by a single scaling factor, or by a value ap- 
propriate to_a_singlecharacteristic wavelength for each reaction in 
question JRawlingsl 19881) . In this study we use a more accurate ap- 
proach that can encompass the variety of situations that we wish to 
investigate: Photoreactions are calculated by fully integrating the 
known cross-sections and oscillator strength s (for line tran s itions) 
over the specified radiation field. Following Ivan Dishoeckl 1 119881) 
we treat the contributions to the photodissociation rate arising from 
continuous absorption and those arising from absorption in discrete 
transitions as separable. These can be written as 






: I a(X)I(X)dli 



and 






Tie 



Kifui^u^tHKi) 



(1) 



(2) 



respectively, where k is the photoreaction rate coefficient, %. is 
the wavelength, ly^ is the flux, Ct(A,) is the cross-section for the 
reaction, /„/ is the oscillator strength and the other symbols 
have their usual meanings. Photoreaction cross-sectional data and 
oscil lator strengths were taken from a v ariety of sources (e.g. 
Ivan Dishoeck 1988: Rawlings. Drew & Ba rlow 1993) as described 



Rollins & Rawlingj <2012l) . We assume that the gas is within 



the carbon-neutral zone (equivalent t o being beh i nd the photodis- 
sociation region: As was shown in iRawlingsl ( Il988l) , complex 
chemistries are unsustainable in ClI regions). We therefore assume 



that H2 and CO are fully self-shielding and that carbon photoion- 
ization is effectively suppressed. The lower limit of the integral 
in equation 1 is therefore set to the ionization threshold for Cl 
(IIO2A). 

This study is too generalised to consider the radiation fields 
of individual sources so we instead opt to characterise the radiation 
field by that of an (appropriately diluted) single temperature black- 
body (Tbb)- In addition, as we have investigated a wide range of Tbb 
which may also be time-dependent, we have utilised a parametric 
fit to the dependence of the photoreaction rates to Tbb- Thus, after 
calculating the photospheric rates for several values of Tbb we fit a 
simple quadratic to the logarithm of the rates: 

log(fc) =a + h. log Tbb + c \og{TBBf (3) 

where the coefficients a.b and c are specific to each reaction and 
are defined for a certain temperature range. The ranges used are as 
specified in Table 4. 

For the geometrical dilution we assume simple spherical sym- 
metry, in which case, at a distance R from the star, the the dilution 
factor is given by 

where Ri, is the photospheric radius of the star. 

To ascertain the viability of the efficient formation of C2H2, 
we will need to study the chemistry in the vicinity of the 'conden- 
sation radius' where dust grains start to form. In the case of AGB 
envelopes, this is also where outflow acceleration due to radiation 
pressure on newly formed dust grains starts. For AGB winds, the ra- 
dius o f condensation, R^, is given bv lCherchneff. Barker & TielensI 
1I1992I) : 



Rc=R.\ 1 



1-2 



4.85 



^ phot J 



-0.5 



(5) 



Substituting this into equation ([4J, the dilution factor becomes 



W=- 1.0- 



1.0-2.0 



-^ 



4.85' 



'p/io/ , 



(6) 



which is only dependent on the ratio of the temperature of the pho- 
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tosphere (Tpi,„,) to that in the gas {Tgas)- In the case of the AGB 
winds, the temperature of the gas parcel is set at lOOOK, which is 
the midpoint of the temperature wi ndow for PAH form ation (900- 
1 100 K) found by Cherchneff. Barker & TielensI ( Il992l) . 

Finally, we include a (small) contribution to the photorate de- 
riving from the loca l internally generated UV field induced by co s- 
mic ray ionizations jPrasad & Tarafdanll983HGredel et al.lll989l) . 



3 PHYSICAL AND CHEMICAL PARAMETERS 

We have investigated four different carbon-rich (C>0) astrophysi- 
cal environments. We have studied acetylene formation in the out- 
flows from: 

(i) Thermally pulsing asymptotic giant branch (AGB) stars, 

(ii) Novae, 

(iii) G/K-type Carbon stars, and 

(iv) R Coronae Borealis-type stars. 

The physical conditions are quite different in these four classes of 
objects, yet they are each known to be prolific sources of carbon 
dust. 

Firstly, we have exten ded the work of 
ICherchneff. Barker & TieTeiisI (Il992h . which modelled PAH 
formation in the post-shock gas in thermally pulsing asymptotic 
giant branch (AGB) stars. These are well known dust production 
sites. Including a more complete photochemistry has allowed 
us to determine more accurately the stability of C2H2 in this 
environment. We have applied our model to examine the viability 
of PAH formation in the ejecta around novae, some of which (but 
not all) are also known to be efficient dust-producers. We have 
also investigated the outflows from stars that are somewhat hotter 
than AGB stars, so that we can ascertain the possibility of soot 
chemsitries operating in these circumstellar environments. Finally, 
we have modelled the conditions required for acetylene-driven 
dust formation in R Coronae Borealis (R Cor Bor) stars, which are 
well-known to be prolific dust-producers. The initial conditions 
selected for each of these four scenarios are discussed below. 

There are considerable variations (and uncertainties) in the ob- 
servationally inferred elemental abundances for each source type. 
However, the C:0 ratio is a key parameter. This determines the ef- 
ficiency of the hydrocarbon chemistry and, hence (whatever reac- 
tion pathway is followed) the viability of carbon dust formation. 
To effect the requisite variotions in the C;0 ratio we have adopted 
cosmic abundance values for all elements other than carbon; 0.85, 
6.0 X 10"^ 4.6 X 10""^, 1.4 X 10"^ 3.5 x 10"^ and 3.2 x 10"^ for 
He, N, O, S, Na and Si respectively. The abundance of carbon is 
then varied to give the (observed) C:0 ratio as specified in Table 4. 



3.1 Case i: AGB Stars 

High carbon-oxygen ratios (> 2) are unlikely in the Galaxy (al- 
though common in the Large Magellanic Cloud) and are ex- 
cluded. L ower ratios of C : O ~ 1.5 are e xpected for most car- 
bon stars JBergeat. Knapik & Rutiivll200lh . iLambert et all l ll988l) 
found slightly subsolar values for the oxygen and nitrogen abun- 
dances in AGB stars, but our approximation of using solar abun- 
dances is essentially valid. This is consisten t with the abundances 
used bv lCherchneff. Barker & TielensI (1 1992h to model an AGB star 
that is similar to the well-kn own source IRC-H10216. We can there- 
fore compare our results to ICherchneff. Barker & TielensI (Il992r) 
and distinguish the effects of including a complete photochemistry. 



In a st udy of of 390 carbon rich stars bv lBergeat. Knapik & Rutilvl 
yOOlj), the highest effective temperature evaluated was 3870 ± 
1080 K. We have therefore investigated photospheric temperatures 
in the range T=1500-3500K. The initial conditions and parameter 
ranges that we consequently adopt for AGB outflows are given in 
Table 4. 



3.2 Case 11 : Novae 

It has long been established that dust is sometimes for med, rapidly 
and efficiently, in nova outflows (e . g. I GallagheJI 1 977r) . In the most 
extreme circumstances, a steep decrease in the visual emission oc- 
curs, accompanied by a strong rise in t hermal infrared emission - as 
an op tically thick dust shell is form ed Jciavton & Wickramasinghd 
Il976l) . Following iRawlingsl l ll988l) . we adopt the following rela- 
tions between the physical parameters: 



Amy = rri\ x fo , 

Lnova = 1-036 x 10^ x my, 

•■phot ~ 9.97 X 10^ X Lnova, 



'phot - 



15280 X W^-rr 



(7) 
(8) 
(9) 

(10) 



where niy is the rate of magnitude decline in the optical, r„jjQj is 
the initial photospheric radius, Lnova is the nova luminosity and Iq 
is the number of days post outburst. 

The values for the rate of visual decline from maximum (ni,. = 
0.03 mag d^' ) and time at sta rt of calculation {tQ= 30 days p ost out- 
burst) were adopted foUowing lRawlings & Williamsl(ll989l) . Novae 
can be classed by the time taken for a nova to diminish by 2 magni- 
tudes below maximum visual brightness. Generally, 'slower' novae 
(i.e. those with smaller values of ra,,) are more effective at produc- 
ing optically thick shells of carbon dust. Parameters for a typical 
dust-producing nova are given in Table 5 and are used in our model. 
The parameter range investigated and initial conditions are given in 
the Table 4. The fractional a bundance of carbon is muc h higher 
than in the case of AGB stars JRawlings & Williamsl 19891) . This is 
due to the thermonuclear processing of the gas and dredge-up of 
white dwarf material that occurs following the thermonuclear run- 
away that drives the nova outbu rst. The upper li mit of density was 



M 



inferred from observations (e.g. iRawlingsll 1988b 



3.3 Case ill: Carbon Stars in Spectral Class G-K 

We investigate stars with photospheric temperatures in the range 
3750-6000K. Low mass main sequence stars are not likely to 
have a high C:0 ratio. However, CH subgiants are low luminos- 
ity peculiar giants, which have a high C:0 ratio in their atmo- 
sphere as a result of mass transfe r from an evolved companion 
l lSmith. Coleman & Lambert 1 19931). The upper l imit of C:0 ob- 
served in these stars is 1.5 jLuck & Bondlll982l) . The parameter 
ranges chosen to represent these objects are given in Table 4. We 
have optimised the parameters for dust production by only con- 
sidering the highest values possible for the C:0 ratio and the gas 
density. 



3.4 Case iv: R Coronae Borealis stars 

R Coronae Borealis stars (RGB's) are yellow supergiants which are 
hydrogen deficient and carbon rich. These stars show irregular de- 
clines in visual brightness caused by the production of thick dust 
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Table 4. Physical and chemical parameter ranges for the four different cases being investigated. 





AGB 


Nova 


Stars of spectral Class G- 


■K 


R Cor Bor stars 


[C/O] 


1.5 


21 


1.5 




7.2 


Density Range (cm"') 


10«-10" 


lO'-lO'l 


lO'-lO" 




lO'-lO'^ 


Photosphere Temperature (K) 


1500-3500 


15,000-25,000 


3750-6000 




4500-8000 


Gas Temperature (K) 


1000 


2000 


1000 




1000 


Dilution Factor 


see eqn (6) 


5.59x10-5 


see eqn J6) 




6.25x10-"* 



Table 5. Table of nova parameters 



Nova Parameter 



30 days post-outburst 



Rate of magnitude decline (optical), my 
Luminosity of nova, Lnova 
Photospheric radius, r_[jQ( 
Photospheric temperature, Tpj^j 



0.03 mag d-' 

31080 

25 Re (=Rp.o) 

20,100 K 



clouds JCrause et al.ll2007l) . There are also variations in observed 
colour that correspond to the mi nimum in visu al brightness and 
which are attributed to pulsations ( lEfimovlll988h . These stars have 
fast stellar winds t hat feed into a thin dust sh ell found hundreds of 
stellar radii away ( IWalkerlll986l : lFeaslill986r) . The optical proper- 
ties of the ejected dust differs from the ISM, with a UV absorption 
peak at 2500 A Ijeffers et al. 2012). Models indicate that this cor- 
responds to small glassy or amor phous carbon grain s (soot) formed 
in a hydrogen-poor environment (IHecht et al.ll984l) . 

The dust formation mechanism is not at all well understood 
and there has been much debate as to where the dust condensa- 
tion occurs. There are two main models: i) the first assumes the 
grains are formed at ~20 R*, and ii) the second proposes that 
dust is formed in the photosph ere of the star and moves quickly 
away due to radiation pressure (IPavne-Gaposch kin 1963). The ra- 
tionale behind the first method is that the stellar surface tempera- 
ture greatly exceeds the gra i n condensation temp erature for carbon 
( lFeastlll986l : lFadevevlll988l) . IWoitke et"aLl ( ll996l) proposed that the 
pulsations stronger than a certain critical amplitude could cause at- 
mospheric shocks to form in a process similar to that proposed for 
AGB outflows. This yields post-shock gas temperatures and den- 
sities of T~1000-1500 K and n~ 10^ - 10^ cm^^ w hich may be 
[ ropit ious for dust formation. In fact, observations by Ijeffers et al.l 
20121) indicate that there are two populations of grains which may 
indicate that both formation processes are occurring. 

We investigate C2H2 formation on the assumption of the first 
model, i.e. molecule/dust formation at a radius of 20 R* from the 
star The te mperature of the gas in th is region is observed to be 
< 1500 K lOoeres & Sedlmavil[l993) . The effective temperature 
of RGB's range from 4.000-20,00 K, with mean values of 6000- 
8000 K for dust-forming sources JLambert & Radii 9941) . The pa- 
rameter ranges adopted for this scenario are given in Table 4. The 
abundances for carbon, nitrog en and helium were obtained from 
iGarcia-Hemandez et al.l 00091) , whereas all other abundances are 
assumed to be solar 



4 RESULTS 

C2H2 formation is a necessary, but not sufficient, criterion for car- 
bon dust formation via PAH chemistry. Therefore we need to estab- 
lish the criterion for a positive result, i.e. the minimum fractional 



abundance of C2H2 that will explain the requisite abundance of 
dust n ucleation sites. To do this, we follow Rawlings & Williarng 
( Il989l) : where the number of grains is equated to the number 
of nucleation sites. Considering the case of a nova, the number 
of grains, Ng, inferred from infra red observations is ~ 7 x 10 
JClavton & Wickramasinghelll976r) . The dust shell condenses at a 
radius re ~ 8 x 10 cm, and has a typical thickness of ~ Jv/lO. 
The shell volume is therefore V ~ 6 x 10 cm , implying a num- 
ber density of grains of ~ 10^^ cm--*. Note that this is a lower 
limit, subject to the ejecta configuration. The gas density is ~ 10 
cm^ , hence the minimum fractional abundance of nucleation sites 
needed is 10^'^. On the assumption that each C2H2 molecule is 
efficiently converted into a nucleation site (i.e. that there are no ki- 
netic bottlenecks in the formation of larger molecules), we can very 
crudely adopt this value as the threshold for the minimum fractional 
abundance of acetylene which can promote a PAH-driven dust for- 
mation pathway. 

Results from our models are shown in Figures [TBI as contour 
plots of the equilibrium abundance of C2H2 as a function of density 
and photspheric temperature. We discuss the results for the four 
cases in sections 94.11- 34.41 



4.1 Case i: AGB 

The results for AGB winds are presented in Figure |T] The figure 
shows that the C2H2 abundance is not strongly sensitive to the den- 
sity. There is, however, a strong dependence on the photospheric 
temperature. For photospheric temperatures below '--^2500K, C2H2 
can be produced with the requisite efficiency for any density within 
the range explored (~ 10 — lO^cm^ ). The lack of senitivity to 
density implies that if C2^2 leads to PAH formation, then the pro- 
duction of dust will not be crtically limited to the geemetrical dilu- 
tion of the wind. 



4.2 Case 11 : Nova 

The results for nova outflows are presented in Figure |2] and show 
the opposite trends what was found for AGB winds; a strong 
dependence on density, but only a weak dependence on photo- 
spheric temperature. Despite the hi gh abundance of carbon in novae 
JStarrfield, Truran & SparkslBOOOl) , the photospheric temperatures 
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Figure 1. Contour plot showing the logarithm of the equilibrium fractional 
abundance of C2H2 in the circumstellar envelope of a carbon-rich AGB star. 
The density and photospheric temperature are treated as free parameters. All 
other parameters are fixed and given values as specified in Table 4. 
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Figure 3. Contour plot showing the logarithm of the equilibrium fractional 
abundance of C2H2 in the circumstellar environment of a star with spectral 
class G-K. The density and photospheric temperature are treated as free 
parameters. All other parameters are fixed and given values as specified in 
Table 4. 
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Figure 2. Contour plot showing the logarithm of the equilibrium fractional 
abundance of C2H2 in the ejecta from a nova. The density and photospheric 
temperature are treated as free parameters. All other parameters are fixed 
and given values as specified in Table 4. 



Figure 4. Contour plot showing the logarithm of the equilibrium fractional 
abundance of C2H2 in the circumstellar environment of a R Coronae Bore- 
alis star. The density and photospheric temperature are treated as free pa- 
rameters. All other parameters are fixed and given values as specified in 
Table 4. 



are ten times higher than for AGB stars and, even after significant 
geometrical dilution is taken into account, C2H2 is highly suscepti- 
ble to photodissociation and never attains appreciable abundances. 
This is a strongly negative result and supports the findings of pre- 
vious studies (e.g. lRawlings & Williamsll 19891) - C2H2 can not be 
produced in irradiated nova environments at abundance levels suffi- 
cient to seed dust formation via PAH chemistry. The only way that 
a positive result could be obtained is if the ejecta were extremely 
dense (n> 10 cm^) allowing three-body chemical reactions to 
take place, or strongly shielded from the nova radiation field. 



4.3 Case ill: Stars in Spectral class G-K 

Results for the circumstellar environment of carbon-rich stars with 
photospheric temperatures of 3750-6000K are presented in Fig- 



ure [3] The critical fractional abundance of C2H2 is not reached. 
However, for the lowest photospheric temperatures (T<4500K) 
and highest densities (n> lO" cm^^) the abundance is very much 
closer to the threshold than was the case for the nova outflows. We 
can conclude that, whilst the spherically symmetric model produces 
a negative result, a positive result may be obtained in the event of 
ejecta clumping and/or stratification. 

However, unlike the situation with AGB outflows, this result 
is not robust to density variations so, if the clumps were to dissi- 
pate within the temperature window for PAH formation then the 
C2H2 abundance (and dust formation efficiency) would be strongly 
diminished. 
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4.4 Case iv: R Coronae Borealis stars 

The results for the gas in a stellar wind located at a radius of 20 R* 
from an R Coronae Borealis star are presented in Figure 

A positive result is obtained for part of the range of the physi- 
cal conditions that we have investigated. The dependencies are sim- 
ilar to those found for the outflows from G-K stars in that there is a 
strong dependence on density and a weaker dependence on photo- 
spheric temperature. The requisite fractional abundance of C2H2 is 
obtained for densities that Js 1 — 5 x lO" cm^-'. For photospheric 
temperatures Tp/„„ < 6000 K, a positive reult can be obtained at 
lower densities. We again note that any density enhancements that 
may exist within the outflow will strongly strengthen this result. 



5 CONCLUSIONS 

We have used a simple model of the steady-state chemistry of sev- 
eral different astrophysical dust forming environments to consider 
the stability of acetylene (C2H2). In previous studies of dust forma- 
tion, based on terrestrial flame chemistries, this molecule is a pre- 
cursor to PAH formation and so its presence at appreciable abun- 
dance levels can be used to estabish the viability of PAH forma- 
tion and, by inference, dust formation based on PAH condensa- 
tion chemistry. However, C2H2 is highly suceptible to photodis- 
soocation by stellar ultraviolet. The novel approach in this study is 
that we have carefully calculated the photodissociation rates for the 
molecular species involved in the C2H2 chemistry, subject to the 
local radiation field. Previous studies had made major simplifying 
assumptions in this respect. 

In each case we have investigated the dependence of the C2H2 
abundance as a function of gas density and photospheric tempera- 
ture. In general, we find that the inclusion of accurate photorates 
significantly restricts the parameter space within which efficient 
C2H2 formation (and hence, presumably, PAH formation) can oc- 
cur. 

In the case of the winds driven by thermally pulsing AGB 
stars, we find that - as in previous studies - C2H2 can be formed 
very efficiently indeed and that a PAH chemistry is therefore at least 
a viable pathway for carbon dust formation in these environments. 
This is an interesting results as some previous studies (which have 
not included photochemistry in thei r modelling) have not been so 
successful at producing C2H2 (e.g. ICherchneff. Barker & TielensI 
ll992l : lHelling et al.l 19961) . However, we find that this positive result 
only holds for photospheric temperatures that are <2500 K. There 
are dust-producing carbon AGB stars wit h photospheric tempera- 
tures ~ 3500K (e.g. lEglitis & Eglitel 19951) . If we were to include a 
more accurate physical model of the outflows, including shocks and 
clumping, a PAH chemistry may also be applicable in these objects. 
Alternatively, other dust formation mechanisms may be operating. 

As compared to the dust-forming environment of AGB stars, 
the ejectae of novae are subject to a much stronger radiation field 
with a significantly higher black-body temperature. The photodis- 
sociation rates for unshielded molecular species, such as C2H2, 
are very much larger than in the othe r sources considered in this 
study (e.g. lRawlings & Williams! 19891) . As a result C2H2 can never 
attain significant abundances. Of course, other factors should be 
taken into account - such as possible variations in the elemen- 
tal abundances in the ejecta - but the extreme clumping required 
is incompatible with the observed spherical coverage of the opti- 
cally thick dust shell in these objects. In this case, it seems that 
other (non PAH-based) formation channels must be operating - 



involving species that can survive the intense radation field (e.g. 
C2, which is partly shielded by the carbon ionization continuum 
iRawlings & Williams! ( 1198^ ). 

CH subgiants have physical characteristics that are closer to 
Main Sequence stars than stars in the Giant phase. They have 
lower luminosities than equivalent sta rs on the Main Sequence and 
are kn own to be chemically peculiar JSmith. Coleman & Lambert! 
Il993!) . having a carbon over-abundance. The winds from these stars 
are not expected to have the same shock-generated density en- 
hancements as in the AGB outflows. Thus, even with the extreme 
parameters we have chosen, we only obtain a marginally positive 
result - with moderate C2H2 abundances obtained at the lowest 
photospheric temperatures and highest gas densities (T<4500K, 
n> 10" cm^ ). Again, an environment that is conducive to C2H2- 
PAH chemistry may exist if significant clumping occurs, but it 
would seem that alternative dust foimation pathways are more 
likely in these environments. 

Finally, we find that C2H2 abundances that are sufficient to 
drive a PAH chemistry can be obtained in the winds of R Coronae 
Borealis stars, provided the density is high enough (^1— 5xlO" 
cm^-'). This sug gests that the large grains detected in the dust shells 
of these objects J leffers et al.l2012n could be produced via a C2H2- 
PAH chemistry operating in dense ejected clouds. By contrast, it is 
possible that the smaller grains observed in these sources may form 
via an alternative mechanism. 

In summary, we find that C2H2 (and PAH) formation is at 
least viable in AGB outflows and the denser parts of R Cor Bor 
winds, but is only possible in the winds of CH subginats if signif- 
icant clumping is present, and is not possible at all in nova winds. 
Therefore, on the basis of these studies, it is evident that the forma- 
tion pathway for (carbon) dust cannot be the same in all environ- 
ments; in particular, the modified flame-PAH chemsitries that pro- 
vide a plausible formation channel in the dust-rich winds of ther- 
mally pulsing AGB stars cannot be applied to environments such 
as novae and may only have a marginal importance in RCor Bor 
winds. It would therefore seem that multiple (and, as yet, largely 
undefined) pathways for the kinetic formation of carbon dust exist. 
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